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An observer, situated several tliousand ligtit-years away from a radio pul- 
sar, finds himself embedded in the diffraction pattern resulting from the 
propagation of the radio waves through the irregular interstellar medium. 
D^' The observer's movement relative to the pattern causes an apparent scintil- 

5-1 . lation of the pulsar. A binary star, situated close to the pulsar's line-of-sight, 

is generating relatively strong gravity waves. 

The rays originating from the pulsar experience a tiny periodic deflection 
/\ • due to the gravity waves produced by the binary star. This deflection dis- 

c^ ■ places the diffraction pattern laterally in a manner familiar from refractive 

interstellar scintillation, except that this gravity wave effect is not dispersive. 
The displacement has the same period as the gravity waves. Its amplitude 
equals the product of the tiny deflection angle and the large distance from 
the binary star to the observer. This periodic displacement can reach a few 
hundred kilometers, which can be comparable to the size of the features in the 
diffractive pattern. Thus, there seems to be a possibility that the exceedingly 
faint gravity waves can manifest themselves macroscopically. 

Observationally, the end effect could be a substantial, deterministic al- 
teration of the scintillation time structure. 



The importance for fundamental physics and for astrophysics of an even- 
tual detection of gravity waves is matched only by the technical difficulty of 
achieving such a detection [1]. 

Most effects through which gravity waves can be detected, such as frac- 
tional proper-distance variations, rates of change in time delays or spatial 
geodesic deflections, have magnitudes of the order of h{D) = H/D, where D 
is the distance from the observer to the source, and iJ is a measure of the 
intrinsic strength of the emission. Thus, the extreme faintness of the waves 
is due as much to the inherent weakness of the emission as to the largeness of 
the distances involved {H << Im and D ~ 10^^ to 10^°m for most galactic 
sources.) 

For example, candidate neutron stars have h{D) ~ 10^^^ in the most 
optimistic estimates {D > lO^^m and H < lO^^m.) Binary stars {D > 
lO^^m and H < lO^^m) could reach in principle h{D) ~ 10^^°, but they 
emit at much lower frequencies. This makes their detection by ground-based 
hardware extremely difficult, due chiefly to seismic noise [2,3]. 

An example of high-frequency, intrinsically strong sources is the merger 
of two compact stars. Here H should reach about 100m, but the expected 
rareness of such events puts most of them at extragalactic distances {D > 
lO^^m.) Still, the resulting amplitudes from the closest events {h{D) ~ 
10—21) are at just about the projected threshold of sensitivity of the LIGO 
detector, which is already at an advanced stage of development [1,2]. 

Recently, an alternative approach to gravity-wave detection was suggested 
[4-7] that tries to circumvent the stringent constraint: effectr^ h{D) = H/D. 
This was attempted by considering effects suffered by electromagnetic rays 
when they propagate passet a source of gravity waves. It was shown that 
if a source of gravity waves lies at a small distance ( "impact parameter" ) b 
from the line-of-sight of an electromagnetic source, the latter will have its 
rays deflected by an angle a^^ ~ 7r^/i(6), and delayed by an amount 6t that 
changes at a rate d6t/dt ~ h{b). 

It was then pointed out that there exist several actual astronomical 
conflgurations where one has b << D, and hence h{b) >> h{D) for the 
same gravity-wave source. A good illustration of this is when a pulsar 
is the electromagnetic source, and a neutron star which is in tight orbit 
around it is the gravity-wave source. For instance, for the binary pulsar PSR 
B1913+16, b comes to only about lO^m every 7^45'. For the optimistic esti- 
mate H fa 10~^m, the expected magnitudes of ray deflection and time-delay 



variation are agu, ~ 10"^" and d6t/dt ~ 10~^^sec/sec. This is not too far 
below today's achievable sensitivities. 

These same waves have an amplitude of only h{D) ~ 10^^^ when they 
eventually reach the Earth {D ^ 1.5 x 10^°m). There is little hope, in the 
foreseeable future, of detecting so weak a gravitational radiation by interfer- 
ing with it at the level of the solar system. Thus, the approach proposed in 
[4-7] amounts to using electromagnetic waves to probe gravity waves in the 
vicinity of their source, where they are still relatively strong. 

Another type of astronomical configuration where the above applies is 
when the source of gravity waves is a binary star that lies close to the line- 
of-sight of the pulsar. For example, the /i-Sco system has H ^ 0.06m and 
D ^ 3 X lO^^m. Hence, the magnitude of these waves, when they reach 
the Earth, is only h{D) ~ 10^^°. If this system were lying at 1" from 
some distant light source, it would induce light deflections and time delays 
comparable to the (very) relatively large values found above for the neutron- 
star case. Furthermore, we shall see, in the following, that the alignment 
constraint is less stringent for the effect described here than for [6,7]. 

In a sense, then, in the cases identified in [6,7] the distance D does not 
suppress the effect of gravity waves {effect ~ h{b) = H/h,) except perhaps 
for the deflection effect, which it can sometimes prevent from translating into 
angular-position shifts through the lens equation. 

The present letter suggests a way of taking this a step further: the very 
largeness of D may contribute to making gravity waves more accessible to 
observation. This could be achieved by exploiting the scattering of pulsar 
radio waves by the interstellar medium [8-13]. 

In what follows, we describe the interstellar scintillation of pulsars in 
heuristic terms that summarize an extensive body of observations and the- 
oretical derivations, about which there now exist some excellent books and 
review articles [14-20]. 

The scattering of pulsar radio waves by inhomogeneities of size a in the 
interstellar electron density result in a "diffraction" pattern for an observer 
in the Fresnel region, i.e., at a distance D > a/6s, where 6s is the scattering 
angle. The scattering is strong or weak according to whether or not the phase 
fluctuation (A0 ^ 2n6sa/X, A being the radio wavelength) is larger than 
about one radian. Hence, the scattering is strong and the fractional intensity 
fluctuations are of order one if the observer is in the Franhofer region: D > 
a^/A. The velocity {V) of the observer relative to the interference pattern 



causes the pulsar intensity to appear to scintillate on a time scale r = S/V , 
where S is the typical size of the features in the pattern. In principle, the 
effect of internal motions in the interstellar medium should be added to that 
of the relative motions of the pulsar and the Earth. 

Observationally, most pulsars show strong scintillation at wavelengths 
A > 10cm. Observations have also allowed turbulent motions within the 
interstellar medium to be ruled out for most (but not all) lines-of-sight of 
interest . 

Over the last thirty years or so, the theory of interstellar scintillation has 
been developed to the point of being able to account for most of the phe- 
nomena displayed by pulsar radio spectra. However, until relatively recently, 
observed long time-scale variations in some spectral characteristics remained 
largely unexplained. Eventually, these slow variations were shown to be also 
due to the action of the interstellar medium, rather than to intrinsic proper- 
ties of pulsars [12,13]. 

Interstellar inhomogeneities with scales a > OgD do not contribute to the 
diffractive scintillation discussed above. Rather, they cause an overall tilt 
of the electromagnetic wavefronts, which are otherwise randomly perturbed 
on the smaller diffractive scale. This phenomenon amounts to a large-scale 
refraction, which consists in ray bundles being deflected by an angle Or- The 
variation of 9r on very large scales (refractive focusing), again combined 
with the observer's relative velocity, causes deep modulations in the pulsar 
intensity over time scales that sometimes exceed several months. 

Closer to the subject of the present paper, refractive scintillation is also 
responsible for the frequency drifting of diffractive features which is clearly 
present in many dynamic spectra. (These are proflles of intensity versus radio 
frequency and time t or spatial position x = Vt.) This is a consequence of the 
dispersive nature of the refraction by the interstellar medium, more precisely, 
of the fact that 6'^ oc A^. 

Now, it is a familiar phenomenon in optics (see e.g. [21]) that changing 
the direction of a bundle of rays that is incident on an irregular medium by 
an amount 6r, shifts the whole pattern laterally by an amount 6sr ~ D6r- 
(We come back to this more explicitly below.) Hence, in the case of refraction 
by the interstellar medium, for which 6r cc X'^, diffractive features at different 
wavelengths develop a relative shift, resulting in the kind of frequency drift 
that has been observed in many dynamic spectra. The formula describing 



the drift rate follows directly from the above by writing 

dt dt ds dOr 2Dd9r 



dX ds dOr dX XV\ 
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where Vj_ is the projection of V unto the direction along which the refractive 
"steering" takes place. 

To the next order of approximation, this steering or group motion of 
diffractive features is accompanied by an alteration (usually slight) of the 
parameters of the diffractive scintillation. A detailed study of these refrac- 
tive alterations can be found in [18]. There, it is shown that the theory of 
refractive scintillation, besides being compatible with most observations, is 
also backed by numerical simulations. 

In the astronomical configurations that we considered earlier (where the 
gravity-wave source is close to the line-of-sight of the pulsar) the gravity 
waves effectively imprint a kink ag^ unto the direction of propagation of 
the narrow ray bundle which eventually reaches the vicinity of the observer. 
Hence, one has a gravity-wave induced refraction with parameters D = Dqe 
and 9r = Ogw Here Dqe is the distance from the gravity- wave source to the 
Earth, and dgyj is the refraction angle corresponding to the deflection angle 
agw. The two are related through the lens equation: 

where Dpc is the distance from the pulsar to the gravity-wave source, and 
DpE the distance from the pulsar to the Earth. 

The initial narrowness of the ray bundle means that, generically, the grav- 
ity waves induce no significant phase gradient amongst the relevant rays. 
(Some possible exceptions will be mentioned towards the end of this letter.) 
Also, the action of the gravity waves decreases rapidly away from the binary 
system. Hence, it is negligible during virtually all the propagation that de- 
termines the diffractive scintillation. Thus, the latter is practically the same 
as in ordinary scintillation. In all, gravity waves generically leave diffractive 
scintillation parameters unchanged (apart from the quasi-rigid shift of the 
diffraction pattern) because the gradient in the refraction angle 6r = 9gu) 
(the amount of focusing) is usually negligible. 

In the context of the present study, the combination of the lens equation 
and the diffraction-pattern displacement equation favors strongly the case 
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where the gravity-wave source is a binary star placed midway between the 
pulsar and the Earth , over the case where it is a neutron star in orbit around 
the pulsar. This is because, in order to maximize the pattern displacement 
OgwDcE, one would like to maximize both Dpc = DpE — Dqe and Dqe- 
Thence, the optimal case is Dqe ~ Dpe/2. 

Before moving on, a word about some of the simplifications made here. 
First, we have tried to keep the discussion to the simplest by considering 
that the observer is confined to a one-dimensional screen, as it is often done 
in the interstellar scintillation literature. Considerations arising from adding 
more dimensions [18] are largely irrelevant in the present context, with the 
exception of the following aspect. Besides being non-dispersive, there is 
another way in which gravity-wave refraction is distinct from its interstellar 
counterpart. When the two polarizations of gravitational radiation are taken 
into account, a given feature in the diffraction pattern will not be oscillating 
along an exactly straight line, but rather along a quasi rectilinear segment 
of a large ellipse centered around the position corresponding to the absence 
of gravity waves. But then, just as in ordinary refractive scintillation, the 
observable part of the effect described here is proportional to the cosine of 
the angle between the (practically rectilinear) shift 6sgu, and the observer's 
velocity V, i.e., it is proportional to V±. 

This is so only because we are considering a ray bundle that is still 
very thin when it propagates passet the gravity-wave source (thickness much 
smaller than one gravitational wavelength) and because then, there is no ap- 
preciable differential of gravity-wave phase, amplitude or polarization over 
a cross section of the bundle. At much larger scales than are relevant here, 
the diffraction pattern is deformed in a deterministic, albeit non-trivial fash- 
ion, due to the extended nature of the bundle cross-section then considered. 
These deformations are closely related to the higher order gravity-wave ef- 
fects mentioned at the end of this paper. 

Finally, we also do not discuss the effect of the non-radiative component of 
the gravitational field generated by the gravity-wave source. That component 
produces a deflection that is much larger than the radiative one, but it follows 
a very different (and precisely predictable) time dependence, as the relative 
motions of the pulsar, the binary star and the Earth steadily change the 
impact parameter b. Hence, it should be possible to subtract that effect, 
along with that of the long time-scale refraction that is due to very large- 
scale interstellar inhomogeneities. 

6 



To see more explicitly how gravity-wave refractive scintillation behaves, 
let us write the slowly varying part of the complex scalar electromagnetic 
field, first in the absence of refraction, as (here we follow closely [18]; see also 
[21]) 

E{x) = ll'^ f dx'Koix - x')e*'^^("') , (1) 



where the function K is the wave propagator: 

Ko{x) =rplexp\i7j:{x/rFoy] . (2) 

rpo is the Fresnel radius rpo = [DpeXY'"^- Iq is the intensity of the incident 
plane wave, x is the spatial coordinate of, say, an intensity peak on the 
one-dimensional screen, (pa is the so-called diffractive phase, generated by 
the random, statistically homogeneous distribution of small-scale interstellar 
fluctuations (a < DpeOs). 

When gravity- wave (or any other) refraction is present, the phase is incre- 
mented by an amount (f)gw{x) which one can represent by a Taylor expansion 
about a reference point Xq on the screen: 

<Pgw{x) = M^o) + 01 (^ - ^o) + <P2{X - Xof + ... (3) 

The intensity Jq is then multiplied by a factor ("gain function") G = 
(1 + Dge^4>2)~^ for DGE^<p2 > —1 (the speed of light is set to one through- 
out.) The pulsar's apparent angular position shifts by an amount Og^ = 
{xgw — x)/Dge, where Xgw is the new position of the peak on the screen, 
i.e. Xgu, = X + 6xguj, where 6xgw = OgwDcE is the quasi- rigid shift mentioned 
earlier. One has 

Xgw = G{x - Xq) +Xo- DgeHi ■ (4) 



The electromagnetic amplitude becomes 

Eg^{x) = {Ghf/^^^ f dx'Kg^iXgw ' X^'^'^^^'^ , (S) 



where now 

Kgnj{x) = Tp^ exp ^m{x/rFf^ , (6) 

with Tp = G^f^rpQ- {^ is an irrelevant phase.) 

For gravity waves, one has generically \Dge^(P2\ << 1, (see below) so 
that G is very nearly equal to one. However, already at the "zeroth-order" 
level (G = 1), a potentially large effect takes place: From eqs.(l-6) and G si 1 
we have, 

EgUx + Og^DcE) ~ E{x) , (7) 

indicating that the pattern has shifted quasi-rigidly along the screen by 
Sxgw ~ OgwDcE, which could possibly measure in kilometers. 

A measure of the observational significance of this effect is the ratio 
6xgw/S, where, as before, S is the typical size of diffractive features in the 
pattern. First, consider an optimal case where a system like /i-Sco, with 
H ^ 0.1m and a gravitational wavelength A^^ = lO^^m, lies at about 1" 
from the line-of-sight of a pulsar. That corresponds to 6 ~ lO^^m if the pul- 
sar is at several kiloparsecs (a typical pulsar distance), and Dqe ~ 10^°m. 
This would imply that ^^p**™-"' r^ lO^^^radian. Then, the diffractive pattern 
could be shifting laterally by as much as (5a;°^**'""' ~ 10.000km. If one relaxes 
the alignment constraint from 1" to 1', the displacement could still be of 
the order of 100km,, a perhaps surprisingly macroscopic manifestation of the 
exceedingly faint gravity waves. 

The spatial scale S, on the other hand, is smaller for longer radio wave- 
lengths, since 

27r^s TeAUeX 

where rg is the classical electron radius and Aug is the electron density fluc- 
tuation for inhomogeneities of size a. 

Observations yield typically 9s ~ 0.2" at A ~ 3m, suggesting that S ~ 
1000km at that wavelength. In practice, A can be increased to about 15m, 
corresponding to S* ~ 200A;m, which is much smaller than 5a;°^**™"'. In all, 
even if the gravity-wave source is either much weaker or much farther from 



the line-of-sight than in the optimal case, the gravity-wave induced periodic 
shift 5xgw could still be of the same order-of-magnitude as the diffraction- 
pattern spatial scale S. 

This concludes our discussion of the lowest-order, and observationally the 
most promising of the effects produced by the gravity waves pointed out in 
[6,7]. 

As pointed out earlier, gravity waves may affect the interstellar scintilla- 
tion at a higher order if it can induce a significant phase gradient within the 
relevant ray bundle. 

Gravity waves alone cannot alter significantly the propagation of electro- 
magnetic waves. Effects such as refractive focusing were shown by several 
authors to be unobservable in the near future [22] . However, when combined 
with the action of the interstellar medium, the effect of gravity waves could 
be made relatively significant in some cases, although not to the extent of 
the lowest-order effect studied above [23]. 

The physics of these higher-order effects is quite involved, and the tech- 
nical details are discussed in a separate publication [23]. Here, we briefly 
outline the motivation for pursuing that direction of research. 

There are two things that a medium has to do to cause scintillation: (1) 
it must bend the rays enough to make them interfere; (2) it must induce a 
substantial phase gradient amongst the rays. Gravity waves are too weak to 
fulfill the first condition by themselves. However, they can, in some cases, 
fulfill the second condition. This can be seen by comparing the modulation 
that they could induce in the optical-path length (AL^^), first with the 
radio wavelength A, and second with the modulation due to the interstellar 
medium (ALj^m-) It can be shown that condition (2) could be satisfied at 
observation wavelengths of a few centimeters. This is due mainly to the 
high-power wavelength dependence of ALj^m, which is proportional to A*^. 

The hope in this picture is, of course, to make the interstellar medium 
fulfill the condition (1), which gravity waves cannot do. This is, however, a 
matter of delicate balance, since at radio waves as short as a few centimeters, 
the scattering power of the interstellar medium itself is considerably weaker 
(recall that Og oc A^.) 

Thus, the question of whether or not the gravity waves, when combined 
with the interstellar matter, can have an observable effect on scintillation 
parameters such as the intensity-correlation frequency bandwidth, the pulse 
broadening or the angular broadening, must be decided by detailed calcula- 
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tion. This involves computing the gravity-wave correction to the two-point 
correlation function of the field (which gives the angular broadening), the 
correction to the two-point frequency correlation function (which gives the 
pulse broadening), and the correction to the four-point correlation function 
(which gives the intensity-correlation frequency bandwidth.) 

One of the difficulties that arise is that one must do the calculations in 
the case of weak interstellar scintillation, as we said. Even in the absence of 
gravity waves, analytical approximations to the four-point correlation func- 
tion were found only in the case of strong scintillation, where it can be 
expressed in terms of lower-order moments. 

What we have just said perhaps sheds some light on why previous re- 
searchers, who looked for observable effects of gravity waves on electromag- 
netic propagation [22], or for gravitational effects on interstellar scintillation 
[24], concluded (correctly) to the negative. In [22], the inhomogeneous inter- 
stellar matter is replaced by (rather than added to) a background of gravity 
waves. Then, as we said, condition (1) is hard to fulfill. In [24], the interstel- 
lar matter is replaced by the stochastic sum of the non-radiative gravitational 
fields of the stars that lie along the line-of-sight of a given pulsar. In this 
case, condition (1) may be fulfilled, but generally not condition (2): The vari- 
ation in the optical path induced by the Schwarzschild metric decreases very 
slowly with the impact parameter [25]: it varies essentially as log(6). This 
is in contrast with the gravity-wave cases of [6,7], where that dependence is 
close to 1/b. This implies that, in [24], the phase differential induced in a 
typical ray bundle tends to be small, hence condition (2) is unlikely to be 
satisfied. 

In conclusion, there seems to be a possibility that the extremely tenuous 
gravity waves can manifest themselves macroscopically. If true, much exper- 
imental work, in terms of noise reduction, would remain to be done before 
such manifestations can be clearly read from the data. In particular, noise 
reduction by exploitation of the high regularity of the waves should be more 
efficient for shorter-period binaries (periods range all the way down to days 
or even hours.) Fortunately, there has been recently a fiurry of activity in 
the field of high- resolution astrometry [26,27], which should greatly enhance 
the chances of finding conveniently located binary stars. Another direction 
of research would be to look for the much more scarce cases where the con- 
veniently located gravity-wave source emits at a much higher frequency than 
considered here. These cases could involve neutron stars or gravitational 
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pulses from violent astrophysical phenomena. In both cases, the frequencies 
could reach lOOHz or more [23]. The effect on the interstellar scintillation 
would then take on a rather different apparence than here, since the velocity 
of the diffraction pattern due to gravity waves could then exceed the velocity 
of the observer. 
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